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We present measurements of AGN type 1 luminosity function in the forbidden line 
[OIII]5007A using data from SDSS DR7. A special attention is paid to NLSyl. A new approach 
in calculating the luminosity function is present. We also account for the large-scale structure 
variations of the Universe density. The results obtained are compared with ones from the liter- 
ature. A prediction of X-ray luminosity function based on our results shows an agreement with 
observations. One of our preliminary conclusions is that NLSyl seems to occupy a more narrow 
range in the nuclear luminosity than BLSy 1, but the average values are within errors. 
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1. Measuring the luminosity function 

At first, we select objects with FWHM(//„) > ^OOkms" 1 (BL AGN) from SDSS DR7 data. 
For our work we use galaxies in the distance interval d c = 55 — 738Mpc. The lower distance limit 
is chosen to avoid the influence of the Local Supercluster [|l|]. The upper limit (z = 0.18) is chosen 
due to impossibility to use our normalization algorithm on larger distances. The main sample is 
divided into two subsamples, NLS (nOOkms^ 1 > FWHM(// a ) < 2000kms~ 1 , 2082 objects) and 
BLS (FWHM(//„) < 2000kms" 1 , 6938 objects). 

In order to evaluate AGN luminosity we used forbidden oxygen line [OIIIJA5007A. It is be- 
lieved that the contribution from star formation to this line is negligible [Q]. Since the emitting 
gas is located in NLR, the AGN luminosity estimates based on this line do not suffer from the 
orientation effects. 

We use a modified V /V nmx method. For each luminosity bin we calculate an average proba- 
bility of detection as a function of comoving distance p(d c ). Knowing this function allows us to 
estimate V max : 



r&c.max 

{Vnmx) = / 4nr 2 p(d c )dr, 

" Ctr.min 



where d c ^ m i n and d c ^ max — lower and upper distance limits of our sample. 

When calculating luminosity function we also took into account variations of density of galax- 
ies per unit of volume due to the large-scale structure. Two facts make it possible: at first, SDSS is 
complete in a wide range of magnitudes. Secondly, the luminosity function of inactive galaxies is 
well known in the Local Universe. 

Moreover, the normalized amount of AGNs in some volume element is: 

N AGN, norm = NagN — " j 

Pgal 

where Nagn is an observed number of AGN, (p ga i) is an average (calculated from luminosity 
function) density of galaxies, p ga i — the density of inactive galaxies which is actually observed. It 
is necessary to consider only galaxies with magnitudes in the completeness interval of the SDSS. 
With help of the luminosity function of inactive galaxies, one obtains: 

{pgal)_L s£h Y_ 

— — ^3 

Pgal i^obs 

where L vc /, is an integral luminosity, calculated using the Schechter function, L () t, s — the actually 
observed total luminosity of galaxies in the considered volume V and luminosity interval. 

The transition from L to p is possible only when the parameters a and L* of the Schechter 
function are fixed. It is the case for galaxies in the Local Universe . 

In fig. [j](7ejfr) we show the result obtained by applying this normalization to a sample of normal 
galaxies split into bins by absolute magnitude. A plateau followed by a rapid decline when visual 
magnitude reaches the completeness limit of SDSS survey is clearly seen. This is exactly the 
expected behaviour, so our normalization is adequate. 



The NLSyl luminosity function 



Andrey Ermash 



6*10 



M r q-19.0,-18.5) 
M r q-19.5,-19.0) 
M r q-20.0,-19.5) 
M r q-20.5,-20.0) 
M r q-2 1.0,-20.5) 
M r q-21. 5,-21.0) 
M r q-22.0,-21.5) 




X 3x10 




6.5 7.0 7.5 

lo 8 L [OIII] /L s U „ 

Figure 1: (LEFT) The result obtained by applying our normalization to a sample of normal galaxies split 
into bins by absolute magnitude. (RIGHT) Luminosity functions calculated in [ODTJ5007A. a) Syl [g] b) 
Sy2 [§] c) Type2 AGN [(| d,e,f) NLSyl,BLSyl,NLSyl+BLSyl - this work. 
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Figure 2: (LEFT) Luminosity functions in B band obtained by converting from [OIII] orH a LFs. Group 1. 
(a,b,c) NLSyl, BLSyl, NLSyl +BLSyl converted from [OIII], this work. Group 2. d) Type2 AGN, [OIII] 
[| e) Syl, [OIII] [| f) Sy2, [OIII] g] g) Syl+Sy2, H a [| Group 3. h) Syl, H a § i) Syl, H a @] Group 
4. g) QSO calculated in the B band [B (RIGHT) a) The observed Syl X-ray luminosity function (0.5-2kev) 
[0. The predicted Syl X-ray luminosity functions, based on b) Our [OIII] LF. c) [OIII] LF from [g]. d) H a 
LF from [§. e) H a LF from @. 



2. The luminosity function 

The obtained luminosity functions of NLS, BLS and Syl (NLS and BLS together) are shown 
in fig. fright). 

For a comparison we also plot the luminosity functions from []2]] and [||] calculated in the same 
spectral line. Note that in [Q] LF is obtained for Syl and Sy2, while in [Q] only for Sy2. 

As can be seen in fig. fright), the results of these two papers are in an agreement with each 
other. Our luminosity function is significantly lower. 

In order to compare our results with as many as possible LFs calculated by different authors 
we convert our LF from [OIII] A 5007 A to the Johnson band B according to [J3|]. 

The result is shown in fig. ^left) where the following LFs are also plotted: 
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1. Converted from [OIII] to M B ([|, g] and this study) 

2. Converted from H a @, |, |] 

3. The quasar luminosity function at redshift 0.4 < z < 0.88 measured in B band [gj. 

As can be seen from these figs LFs obtained by different authors form two distinct groups 
(groups II and III on fig. ^left)). The first and the second ones consist of the results of [g, [J and 
[Q, |5p, respectfully. LF of Sy 1 obtained by us is between these contradicting results. It is important 
that our LF agrees well with the LF of quasars [||] on highest luminosities. 

It is remarkable that <p(NLS)/<p(BLS) is not constant. Instead it is actually a function of the 
luminosity (fig. fright)., bottom panel). It has a maximum between log (Lr o/// i/L ) = 6.8 and 7.6. 
At larger luminosities it slightly decreases. The decrease towards lower luminosities is ~ 0.5dex to 
log (Loiii/Lq) = 5.5. 

A putative explanation can be proposed if the results of [Q] are considered. Authors studied 
samples of NLS and BLS selected from SDSS DR3. Ones of the key results of their work are the 
distributions of log (L/L e dd) and log (Mbh/Mq). The normality of this distributions implies that 
log (Lboi) also has a normal distribution with the following parameters: 

• NLS:jU = 44.589 ± 0.052, a = 0.404 ± 0.043 

• BLS:jU= 44.535 ±0.088, a = 0.560 ±0.061 

The average values are equal within the errors. But the widths of the distributions differ 
significantly. Thus, <p (NLS)/<p (BLS) ratio has a maximum at log (L& /) = 44.64, which corresponds 
to log {L\ohi\/Lq) = 7.4. The dispersion of NLS distribution is lower, thus NLS occupy a narrower 
range in luminosity than BLS. 

3. The x-ray luminosity function 

We calculate a predicted soft x-ray(0.5-2kev) luminosity function based on the L\ m\ data. 
[OIII]A5007A luminosities are converted into the soft x-ray following H^. The result is plotted 
in fig. fright). For a comparison, predictions based on the H a luminosity functions from [Q, [l^] 
and [OIII] from |Q] are also plotted in fig. fright). H a luminosities are converted to the bolometric 
ones following [Q], then from the bolometric ones to Lo.5-2ifcev following [||]. The x-ray LF of 
Syl based on the ROSAT, XMM-Newton and CHANDRA data [JTTJ] is also shown in fig. fright). 



It should be stressed that our prediction corresponds to the actually observed x-ray luminosity 
function better than LFs from other studies. 

4. Results 

• A method of evaluation of AGN luminosity function based on emission-line data accounting 
for variations of the density of galaxies due to the large-scale structure 

• The [OIII]A5007Aemission-line luminosity functions for NLS, BLS, Syl (NLS+BLS) are 
obtained. 
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• At largest luminosities the LF of Syl agrees very well with the one of quasars [Q|. 

• A prediction of soft x-ray luminosity function based on our emission-line LF is in an agree- 
ment with the observed x-ray LF 

• (j)(NLS)/(j)(BLS) is not constant being a function of luminosity. It has a maximum between 
log(Lr o/7/ i/L ) = 6.8 and 7.6. This agrees with results of [|7|]. The average luminosities of 
NLS and BLS are equal within errors, but widths of the distributions differ. NLS occupy a 
narrower range of the AGN luminosity. 
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